Observations of cluster substructure using weakly lensed sextupole 

moments 
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Since dark matter clusters and groups may have substructure, we have examined the sextupole 
content of Hubble images looking for a curvature signature in background galaxies that would arise 
from galaxy-galaxy lensing. We describe techniques for extracting and analyzing sextupole and 
higher weakly lensed moments. Indications of substructure, via spatial clumping of curved back- 
ground galaxies, were observed in the image of CL0024 and then surprisingly in both Hubble deep 
fields. We estimate the dark cluster masses in the deep field. Alternatives to a lensing hypothesis 
appear improbable, but better statistics will be required to exclude them conclusively. Observation 
of sextupole moments would then provide a means to measure dark matter structure on smaller 
length scales than heretofore. 



PACS numbers: PACS: 98.65-r, 95.35.-|-d 



I. INTRODUCTION 

The observations of the Supernova Cosmology Project 
and the High-z Supernova Search team as well as 

the Cosmic Microwave Background observation 0| , later 
confirmed by the Wilkinson Microwave Anisotropy Probe 
(WMAP) j5|], suggest the possibility of an accelerating 
expansion of the Universe, which could mean the pres- 
ence of a positive cosmological constant or dark energy. 
These observations give rise to questions as to the na- 
ture of this dark energy /cosmological constant, possible 
predictions for the evolution of the Universe based on 
this assumption, and possible ways to improve or find 
new observational methods that could lead to a better 
understanding of this evolution. 

Weak gravitational lensing methods (for review see 
[1,13,13 and references therein), that allow one to inves- 
tigate the evolution of matter clustering and the growth 
of large-scale structure 0, 0|, are a way to probe 
both dark energy and dark matter. In a sense it is 
a unique way to investigate both the past and future 
of the universe. The large structure growth is defined 
by the growth of the matter density fluctuations pre- 
dicted by inflationary cosmology and dark energy evolu- 
tion (for review see The structure of density fluctu- 
atio ns p rovides information about inflationary scenarios 
[l^ Il4l . On the other hand, the matter distribution 
measurements give bounds on the dark energy equation 
of state, allowing one to predict the future of the universe 

mill Hi. 

The traditional weak gravitational lensing techniques 
0, despite observational superiority in measuring the 
large clumps of matter such as clusters of galaxies (visible 
or dark) with 10^^ are not sensitive to the substruc- 
ture of such clusters or smaller groups and clumps of 
matter. We have investigated a possibility to use more 
sensitive methods to expose such substructure and de- 
tect the presence of smaller clumps. The development of 
these methods and their application to analysis of data 
from the future observational projects like the SuperNova 
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/Acceleration Probe (SNAP) [Tl,|2a, and the Large Syn- 
optic Survey Telescope (LSST) j22| promise to expand 
our knowledge of large scale structure to scales unreach- 
able by other techniques such as microwave background 
anisotropy or traditional weak lensing. These observa- 
tions could make an important contribution to the un- 
derstanding of dark matter structure as well as possible 
inflationary scenarios. 

Arclets are a familiar strong- lensing phenomena |^ 0] ■ 
They are an example of the general property that the 
nonlinear 1 /r deflections of a light stream passing a mass 
concentration will produce, relative to the stream cen- 
troid, a full complement of moments. The curving seen 
in the arclet can be understood as the correlated super- 
position of a quadrupole and sextupole moment, with the 
length of the arclet usually determined by the strength 
of the octupole moment. 

To date, weak lensing has concentrated exclusively on 
quadrupole moments - ellipticity 0, because usually all 
other lensing-induced higher moments are smaller than 
the lensing-induced quadrupole moments (that are al- 
ready small compared to background galaxy ellipticities) . 
A typical measurable cluster has a mass of 10"'^'* solar 
masses and a radius of 500 kpc. Since the strength of the 
quadrupole kick is proportional to the mass and falls off 
like 1/r^, one could get the same quadrupole moment in 
a light stream positioned 5 kpc from an object of 10^° 
solar masses. In the latter example, since the sextupole 
moment varies as l/r"^, the sextupole moment becomes 
100 times stronger, becoming larger than the intrinsic 
sextupole moments of background galaxies. It occurred 
to us that if clusters or groups contain an abundance 
of lower mass clumps, the light streams passing through 
these groups might occasionally pass close to a small mass 
clump producing an observable sextupole moment in the 
image Q^. In this case one can expect a correlation 
between the direction of the quadrupole and sextupole 
moments. 

To investigate this hypothesis we considered the im- 
ages of the North and South Hubble Deep Fields (HDF) 



2 



|23l l24j as well as some other Hubble images of known 
clusters, such as CL0024. We used SExtractor software 
[25I I to select either faint images (typically 23<m<29 ^) or 
used z-catalogs (when available) to identify distant galax- 
ies with z>0.7. The "curved" galaxies we sought were 



FIG. 1: Hubble image of galaxy cluster CL0024. 



grees of a quadrupole minimum. We will refer to galaxies 
which have a quadrupole maximum aligned with a sex- 
tupole maximum as " aligned " galaxies. 

Figures Hand 121 show the CL0024 cluster and the loca- 
tion of the background galaxy sample and its "curved" 
members. The arrows point in the direction of the scat- 
tering center and their length is proportional to the 
strength of the quadrupole moments. There is clear ev- 
idence of clumping toward the center of this cluster. To 
our surprise, "curved" galaxies were also clumped in both 
the north and south fields of the HDF survey. The prob- 
ability for the observed clumping to occur by chance was 
about 1% for each field. 

We have carried out other tests that support the hy- 
pothesis that the observed clumping comes from lens- 
ing: 1) "aligned" galaxies were predicted and found to 
be clumped as strongly as curved galaxies, 2) galaxies 
halfway between "aligned" and "curved" were predicted 
and found to have no clumping, 3) "curved" galaxies were 
predicted and found to have smaller moments than all 
galaxies, and 4) "aligned" galaxies were predicted ^ and 
found to have moment strengths distributed as all other 
background galaxies. 
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FIG. 2: m-selected background galaxies of the cluster CL0024 
are indicated by blue dots and foreground galaxies by yellow 
circles. The small stars with arrows indicate the "curved" 
background galaxies with the arrow pointing toward the scat- 
tering center and its length proportional to the strength of 
the quadrupole moment. 

identified as those whose sextupole moment was oriented 
so that one of its minima was aligned within a few de- 



We have considered causes other than lensing for the 
spatial clumping of curved galaxies, such as instrumen- 
tal effects, computational effects, or other physical phe- 
nomena. For example, since the background galaxies are 
known to be clumped, our result could be the conse- 
quence of the fact that, for some unknown reason, galax- 
ies in some groups tend to be more curved than in other 
groups. Or perhaps galaxies of a certain epoch tend to 
be more curved. The tests we have constructed and their 
implications for each alternative hypothesis are discussed 
and the results summarized in Table I. 

Finally, assuming we are indeed seeing lensing, we have 
made an initial attempt to deduce the mass of the objects 
doing the lensing and the mass of the group in which they 
reside. The mass of the group is easier to estimate than 
the mass of the lensing objects themselves, in that the 
former depends only on the i) the size of the observed 
group, ii) the fraction of background galaxies that are 
lensed, and iii) the average induced quadrupole moment. 
The results we present here are preliminary. 

Since in some cases the footprint of the observed galaxy 
could pass through a dark matter clump, we present an 
approach to this general situation as well. While indi- 
vidual mass measurements remain out of reach, it could 
be expected that modeling would supply a basis to de- 
duce well-defined statistical results when analyzing larger 
pieces of sky. 



^ Our threshold setting, typically 10 times the mis of sky noise 
floor, causes objects to appear dimmer than the total integrated 
luminosity magnitude. 



^ We thank B.J. Bjorkcn for pointing this out to us. 
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II. 



THEORETICAL FRAMEWORK 



FigElshows rays from a distant background galaxy be- 
ing deflected by a mass concentration. The apparent im- 
age, as seen by the telescope, is defined by an intensity 
function which depends only on the angle of each ray 
as it enters the telescope. Following a ray backwards, 
toward the apparent image, it is deflected by mass dis- 
tributions, but is known to depart somewhere from the 
source galaxy. It will have a definite position and an- 
gle at the source galaxy (measured relative to the po- 
sition and angle of the centroid ray). This map from 
telescope variables to source variables will be symplec- 
tic since the light geodesies are described by a Hamilto- 
nian. From the point of view of the galaxy, all rays traced 
backward from the telescope come from the same point. 
Hence the two angles at the telescope uniquely describe 
the trajectory through space and the initial position xs 
and ys (and angle x'g and yg) at the galaxy. So the 
"backwards" map can be written as a set of 4 functions 
xs{xT,yT), ysixT^yr), x's{xT,yT), and y'sixT^yr), 
where "T" designates "telescope" and "S" designates 
"source" . The coordinate system for both the source and 
telescope images can be taken to be the pixel grid on the 
focal plane. The position of the centroid trajectory is 
taken to be the origin for both images, i.e. the dipole 
kick suffered by the image as a whole is ignored. Only 
under special circumstances, such as Einstein rings or 
point-to-point focusing, will a ray leaving the telescope at 
two different angles arrive at the same point on the source 
galaxy. We will not need to consider such cases and hence 
will be able to drop the two functions x'g and y'g. We will 
be able to assume that the determinant 



dxs dxs 

dxT duT 

dys dys 

dxT dyT 
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(1) 



The two functions xsixr^yr) and ysixT^yr) can be 
combined into one complex function by defining ws — 
Xs + iys- This complex function can be written in terms 
of the variables wt = xt + iyr and wt — xt — iyr by 
substituting xt — \ {wt + wt) and yT = Yi ~ ''^^)- 
The map equations can then be written as the single func- 
tion ws{wtWt)- Since the transverse width of the light 
stream will be small compared to characteristic dimen- 
sions of the variations of the mass distributions, we may 
expand this function in a power series about the stream 
centroid: 



wsiwT, Wt) — Wt + 



n.m— 



OinmWTWT ■ 



(2) 



The 1 -|- oio combination of terms is a rotation and 
scaling, the aoi term is a quadrupolar distortion, the ao2 
term is a sextupolar distortion, the aos term is an octupo- 
lar distortion, the 020 term is a cardioid-like distortion, 
and the an term is an r^-dependent translation of cir- 
cles, and so on. We will be concerned with the terms 



floi, ao2, ao3, and 020 and refer to them more simply by 
the letters, a, 6, c, and d respectively. We have intro- 
duced the tilde symbol to alert the reader to the 
fact that these coefficients have dimensions, and to dis- 
tinguish them from dimensionless partners we will intro- 
duce later. For a map arising from a single kick agi is 
necessarily real and an = 2a2o = 2d. 

Since the largest mass distribution size (~500 kpc) 
is still small compared to typical path lengths ('^1000 
Mpc), and since the light-deflection angles are small (< 
10~^ radians) and can be calculated by multiplying the 
deflection angles of a non-relativistic particle by 2, one 
can integrate the transverse component of the 1 /r^ force 
from a point distribution along a straight path to get the 
deflection angle 



Ar' = - 



AMG 



AMG 
ro 



-1 



dr 
ro 



AMG 
6r 



1 



(3) 



where Sr = r — tq. 

The potential for this kick is 2^s = AMGLn[r], 
which is the Green's function for the 2 dimensional 
Laplace equation, V'^<f> — AnGp, where p is the 2D 
(longitudinally-integrated) density function for the mass 
distribution. 

The usefulness of the complex variables originates in 
part from the fact that the solution to Laplace's equa- 
tion in empty space can be written as the real part of 
an analytic function. For the point mass source (which 
would be the same as outside a symmetrical distribution) 
one may write 2^s — AMG K ( Ln [x + iy])- Expanding 
about the centroid ray located at some wq = xq + iyo 
results in (a constant is dropped) 



2$A 



-AMG 3? 



Wo 



(4) 



It is useful to introduce derivative operators 



a 



dy 



and -M-r = 

aw 



dx 



dy 



which have the 



property that |h ^ |h ^ 1 and |S = |^ = 0. In 

f f J Qw aw aw aw 

terms of these operators the kick from the potential can 
be written as (the constant dipole kick is dropped) 



l^w's = -2 



3(2$^) _ AMG 
dw wq 



^2V ^0 



(5) 



The first three terms of this expression are the 
quadrupole, sextupole and octupole, respectively (see 
Fig^. The geometry of FigEl shows that Aws = 
Dls Aw'{w,w), where D^s is the angular-diameter dis- 
tance between the source and the lens. And if one wishes 
to obtain the map as a function of the telescope vari- 
ables rather than the lensing plane variables one may 
substitute w = -S^wt- 
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FIG. 3: A diagram showing 3 light rays from a source image 
scattered by a single concentrated mass. The relationship 
Dts^Gts = Dls^Sls expresses the observed displacement 
angle in terms of the deflected angle. Note that the image in 
the radial direction is narrower than the source width. 

Deflection Angle 

-0.1 




Quadrupole Sextupole Octupole 





FIG. 4: The curved line in this graph shows the 1/r kick from 
a point source. At a point ro this curve may be expanded 
in a power series (see eq. |2J yielding linear, quadratic and 
cubic terms that are shown by a tangent line, a dashed curve 
and a dot-dash curve, respectively. The deformation of a con- 
stant intensity circle in a light beam centered at ro by each 
of these terms is also shown here. The straight lines indicate 
the radial direction from the origin to the center of the light 
beam. The orientation of these images is such that there is al- 
ways a deflection minimum toward the scattering center. The 
quadrupole plus sextupole yields a slight "banana" shape with 
the radius of curvature originating at the scattering center. 



It follows that the coefficients a, b, and c of this point- 
source map are given by 

^ Dtl ^MG 

as = ^LS-f:. To—, 

Dts Wq 



bs 



cs 



n fDTLyiMG 
Uls 7^ — —r?-. 



\Dts J 



UlS 7^ —=!-■ 



\Dts J 



(6) 



The magnitudes of interest will be the rms value of the 
dimensionless quantities as, b — bga, and c — cs<J^, 
where a is the rms size of the image. Each of these co- 
efficients decreases in magnitude from the previous by 
^g/i^O: where re = j^^^ a is the footprint radius of the 
background galaxy at the lensing mass plane and rg is the 
distance from the center of the footprint to the center of 
the lensing mass. 

For a potential derived from a solution of Laplace's 
equation for a general mass distribution, one can expand 
the potential about the centroid ray to get its local power 
series expression. The map obtained from the derivative 
of this power series potential will necessarily be of the 
form of equation Since the potential is a real func- 
tion, its expansion in terms of w and w will have some 
constraints: 

1. the coefficients of the linear w and w terms must 
be the complex conjugate of one another, so there 
is one complex parameter, which is the magnitude 
and direction of the dipole kick; 

2. the coefficients of the w^and w'^ must be the com- 
plex conjugate of one another, and represent the 
magnitude and orientation of the quadrupole kick; 

3. the coefficient of ww must be real (the coefficient 
of this term must be proportional to -^-^^\q = 
J V'^'&Iq = ttG p\q and hence is zero except when 
the light path passes through a distribution ) and 
represents a magnification; 

4. the coefficients of the w'^ and must be the com- 
plex conjugate of one another, and represent the 
magnitude and orientation of the sextupole kick, 

5. the coefficients of the w'^w and ww'^ must be the 
complex conjugate of one another, and represent 
the magnitude and orientation a third order term 
which is also necessarily zero except when the light 
path passes through a distribution. In this last 
case, the kick can be written in the form Aw — 
dw^ + 2dwTWT-^ d will be proportional to the 



The complex notation is awkward in this case. For real d the po- 
tential function is 3> ~ —dr"^ cos 9, from which Ar' ~ +3dr^ cos 9 
and rA9' ~ —dr'^sm9. 
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first derivatives of p, in fact 



whence 



LS 



D 



TL 



D 



TS 



dw 



6 « - 



AM; 



30 



(12) 



The final element we will need in a minimum the- 
oretical framework is a method to deduce the coef- 
ficients of the map from the image. This process 
begins by noting the surface brightness relationship 
is{xs,ys)dxsdys = irixTyUT) dxrdyT, expressing the 
fact that if the area element is transformed according 
to the map, the number of photons leaving the source 
in that area will be the number observed in the image. 
Defining the (unknown) moments of the source through 
^nm = I Wsw"^ isixs,ys)dxsdys, and transforming to 
telescope variables yields 



nm 



(7) 



Wg{wT, wt)w'^{wt, WT)iT{xT, yT)dxTdy' 



'T- 



This equation can be used to determine an expression for 
the quadrupole map coefficient: 



Mio 



J {wt + awT)^ irixT, yT)dxTdyj 



(8) 



Using equation ^ for a and introducing the rms radius 
of the galaxy to obtain a dimensionless quantity we have 



ba 



-AM: 



3 AM 



30 + 77 



20 



2 Mfi 



3MJ, ' 



(13) 



where a 



i/M/]^. Here we note the interesting fact 
that even with 6=0 the quadrupole term in the map 
can induce a change in the sextupole moment AM30 if 
Mji ^ 0. 

An M21 could originate in the background galaxy or 
be generated by the 3'"'' order kick proportional to the 
derivative of p. To examine the latter possibility we look 
at 



M^ 



{wt + Aw) {wt + Aw) ixixT, yT)dxTdyj 



with Au) = awT 
the equation 



Mfi = Mji 



2dwTWT 



dw^. 



2aM. 



21 



aM. 



30 



5dAL 
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The result is 



(14) 



The moments on the right hand side of this equation 
can be deduced from the telescope image. One obtains 
a quadratic equation for a: one can assume = 0, 
or one can find a statistical ensemble for a in terms 
of an assumed statistical ensemble for M^. Note that 
wH} = so Af]^ is the mean square radius of the tele- 
scope image. We normalize surface brightness functions 
so that J irixT, yr) dxrdyr = 1- 
Solving for a, 



AM 



20 



2Ml, 



2MT^ ' 



(9) 



where we introduced the notation: 
AM„„ = Ml 



M: 



Since [M^q/2M^j] « 0.15, the last term in ® is typ- 
ically only a 2% correction. Though equation ijHl may be 
solved exactly to get 



AM- 



20 



2M/i 



AAf2o '^l 



(10) 



we will be content to approximate the 2'"* factor by unity. 

To obtain an equation for the sextupole moment con- 
sider 



/( 



Wt + awT + 



iTixT,yT)dxTdyT 
3a Mji -\- ... 



(11) 



d — da 



(15) 



-AM21 



AM21 
Ml, 



AM- 



20 



2Ml 



Mi 
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d has the direction opposite to the density gradient and 
in typical cases would be expected to be along the axis 
of the minima of the induced quadrupole moment. 

Finally, since the telescope image has been blurred by 
the point-spread function, one needs to know the mo- 
ments of the point-spread function, and if possible correct 
its systematic effects on the properties of the measured 
image. If the complete (no thresholding) image was avail- 
able one could compute (""" indicates convolved with the 



PSF, and Axt = xt 



etc.) 



M'j.„, = I w^Wt '^TixT,yT)dxTdyT (16) 
w^w^ pirr — t't) iT(frp)dxTdyTdx'Tdy'rp 

{AwT + w't)^ (Awt + w't)™ X 
p{ArT)iT{^T)dAxTdAyTdx'rpdy'rp. 

When the binomial powers are expanded, the double in- 
tegral reduces to the sum of a product of single integrals. 
One finds, for example. 



MJq — M20 



2MfoMf„ 



Mio 



Mio + Mio- 



(17) 
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The left hand side can be determined directly. The first 
term on the right hand side is the quadrupole moment 
of the point-spread function, which can be determined 
presumably by looking at star images. The dipole terms 
in the central product are both zero by definition of the 
centroid. The last term is the sought after quadrupole 
moment of the pure image. Note that if the point- 
spread function has no quadrupole term the blurred im- 
age has the same quadrupole moment as the original im- 
age. However the map coefficient a also involves the mo- 
ment MH, which is different: 
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(18) 



This is a statement that the rms of the final image is 
equal to the rms of the original image plus the rms of the 
point-spread function. For completeness, for remaining 
moments of interest we have 




^4 = ^^3^0 



Ml 



Mi 



Ml, 



(19) 



(20) 



FIG. 5: Two galaxy images (contour plots and 3-D plots of 
surface brightness from the Hubble north field). 



and 



Ml + 2MlMl + Ml + MlMl + mImI 



Ml + 2MfiMi\ -I- M. 



(21) 



In all cases, if the moments of the point-spread-function 
are known, the moments of the original image can be 
found from the smeared image. However this result does 
not strictly hold when the image is taken to be only those 
pixels whose count number exceeds some threshold. We 
will not discuss this problem further here. 



III. GALAXY SELECTION 

The software SExtractor was used to select galaxies 
from the Hubble deep field and to specify which pixels to 
include in the image. Galaxy images were transferred to 
the Mathematica programming environment (see Fig. 
Galaxies with more than one maxima were eliminated. 

Two SExtractor input parameters proved to be signif- 
icant: i) the threshold setting, which we typically took 
to be 10 times the noise floor, and ii) the convolution 
matrix for pixel selection, which we typically took to be 
3x3 or a delta function. If the convolution matrix was 
too broad (larger than 3x3 pixels), important pixels were 
dropped resulting in a change in the computed moments. 
That the dropped pixels were important was deduced 
from the fact that the clumping result described below 
becoms less significant. 

The high threshold is important to observe clumping. 
In fact the clumping result become less apparent below 
5 a. Visual examination of galaxy images reveals that 



at lower thresholds the edges are indeed less crisp, and 
because the sextupole moment is weighted by r"^, this 
moment can be dominated by edge effects. 

Setting the threshold at 10 a noise floor results in im- 
ages which are galaxy cores. About half of all faint galax- 
ies are lost in this cut. Since some fainter galaxies also 
have "crisp" cores, there may be ways to recover those 
by imposing a threshold that scales with core brightness. 




0. 51. 1. 52. 2. 53. 3. 54. 4. 55. 5. 5G. 6. 57. 7. 58. 8. 59. 9. 51 



FIG. 6: Distribution of galaxy rms core sizes. 



Our galaxy-core images have a small rms radius of 
about 2 pixels. The distribution of rms radii is shown 
in Fig. El The distance across the image is typically 5 
pixels. For comparison the rms PSF radius for the Hub- 
ble deep fleld is 1.7 pixels. However for high thresholds 
the point-spread-function spreads the image less than in- 
dicated by its rms radius. 
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IV. 



MAP COEFFICIENT STRENGTHS 



We have carried out studies for galaxies that have been 
selected by magnitude or by z- value ^. Results for all 
quantities were similar. Fig. [71 shows the measured dis- 
tribution of the quadrupole a coefficient for magnitude- 
selected north field galaxies with 



ao 



2MS ' 



(22) 



{M2Q is assumed to be zero, indicates smeared quanti- 
ties are used). The x and y components of the coefficient 
are well fit by a Gaussian distribution with aao= 0.14. 
The distribution for the magnitude of a is the product of 
two Gaussians integrated over angle (a Rayleigh distri- 
bution) characterized by the same &aO. The quadrupole 
moments of selected galaxies were also calculated using 
a "weighting-function" technique. The two methods give 
similar results for the a coefficient. ^ Fig. |S1 shows the 




0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4 0.44 



ao 



FIG. 7: Distribution of the dimensionless quadrupole mo- 
ment coefficient, laol. 
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* We used the z-catalogs from www.ess.sunysb.edu/astro/hdf.html 

and bat.phys.unsw.edu.au/ fsoto/hdfcat.html. 
^ We thank David Wittman for providing this result. 



FIG. 8: Distribution of the dimensionless sextupole moment 
coefficient, |6o|- 

measured distribution of the dimensionless sextupole co- 
efficient 



bo 



Ml 



ML 



(23) 



{M^Q is assumed to be zero, the AAf^ term is ignored) 
for m-selected galaxies in the north field. The x and 
y component of the coefficient is again Gaussian with 
(Tho=0.034 being a factor of 4 smaller than the quadrupole 
moment cTao- Fig. shows the measured distribution for 
the sextupole-order d coefficient 



do = 



ML 



(24) 



( is assumed to be zero and the AAI^q and AM20 
terms are ignored ) for m-selected galaxies in the north 
field. The x and y component of the coefficient are again 
Gaussian. tTdo=0.009, is a factor of more than 15 smaller 
than the quadrupole moment aao and about a factor of 
4 smaller than (Tf,o- Referring to equation H13|l we see 
that the magnitude of the 2'"^ correction term can be 
estimated to be typically a factor of 10 smaller than the 
leading term. In equation (|15|) we can estimate that the 
2"'' and 3'"'' term are typically a factor of 3 smaller than 
the leading term. 
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FIG. 9; Distribution of the dimensionless coefficient Idi 







Finally, Fig. ^| shows the distribution of the dimen- 
sionless octupole coefficient 



Co 



4Ml 



(25) 



It is not much smaller than the sextupole moment: aco— 
0.026. We note that a bi-Gaussian intensity distribution 
with unequal major and minor axes will have a significant 
octupole moment aligned with its quadrupole moment. 
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FIG. 10: Distribution of the dimensionless octupole moment 
coefficient, |co|. 



FIG. 13: The distribution of quadrupole-octupole orientation 
in the HDF North for 324 z-selected galaxies). 



V. MAP COEFFICIENT ORIENTATION 
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FIG. 11: Distribution of the smallest angle btween a 
quadrupole minimum and a sextupole minimum for 324 galax- 
ies in the north HDF. 

Correlated Moments Anti-correlated Moments 




In addition to the magnitude of the sextupole moment, 
its orientation with respect to the quadrupole moment is 
of particular interest to us. A plot of the relative ori- 
entation of these moments is shown in Fig. ^2 for 324 
z-selected galaxies from the Hubble north field. An an- 
gle of 0° indicates that one of the sextupole minima is 
aligned with a quadrupole minimum. We will call such 
galaxies "curved" . When the angle is 30° one of the sex- 
tupole maxima is aligned with a quadrupole maximum. 
We shall call such galaxies "aligned" see Fig. ^1 . 

The orientation of the octupole moment with respect 
to the quadrupole is shown for the same sample in Fig. 
1131 The orientation that results from an induced kick is 
at 0°; the orientation that occurs naturally, and would 
be present for example in a bi-Gaussian distribution is 
oriented at 45°. Though the octupole story is of some 
interest, and the small bump at 0° in the angular distri- 
bution is tantalizing, statistics at this time are too small 
to draw meaningful conclusions and we will not discuss 
the octupole further. Fig. [TUshows the orientation of the 




FIG. 12: Sketches showing the sum of quadrupole and sex- 
tupole moments when moment minima are aligned ( "curved" ) 
and when moment maxima are aligned ("aligned"). 



ita. 
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FIG. 14: The distribution of the orientation angle for do (left 
plot) and d (right plot). 0° is defined by the direction to 
the scattering center as determined by the quadrupole and 
sextupole moment. 



d coefficient with respect to the quadrupole coefficient. 
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VI. CLUMPING OF CURVED AND ALIGNED 
GALAXIES 

The interesting observation about "curved" galaxies is 
that they seem to be clumped. There are two ways to 
study this. The first method is to calculate the num- 
ber of curved neighbors of each curved galaxy that lie 
within a certain distance and compare this to a distribu- 
tion of the same number {Nc) of randomly chosen galax- 
ies. We draw a circle of a fixed radius about each curved 
galaxy in the field. If the circle intersects the bound- 
ary of the field we drop the galaxy from consideration. 
Otherwise we count the number of neighboring galaxies 
within each circle, and make a histogram showing the 
number of galaxies with one, two, three, four, and so on, 
neighboring galaxies inside the circle. To judge whether 
the observed distribution is unusual, we have compared 
this histogram with the average histogram from 300 sets 
of Nc randomly chosen galaxies. Considering that the 
z-distribution of galaxies might play a role in the spatial 
clumping, we decided to require the randomly chosen 
samples to have approximately the same z-distribution 
as the curved set of galaxies. An example of such a his- 
togram is shown in Fig. ^|for the North HD field. We 
call a galaxy "curved" if the angle between quadrupole 
and sextupole minima is between 0° and 8° (see Fig. [TTl 
and Fig. I12|l . Fig. 1161 shows four such histograms for 



R = 320 pixels 



R = 340 pixels 




FIG. 15: A "neighbors" histogram for curved galaxies for 
HDF North in the circle with r = 360 pixels. The background 
(gray) bins represent the averaged random set, the foreground 
histogram corresponds to a curved set. 

several circle radii (in number of pixels). In this radius 
range (320 pixels to 380 pixels) one sees that the curved 
galaxy set consistently tends to have larger numbers of 
neighbors within its circles.^ 



The Hubble deep field images have a drizzled pixel size of 0.04 
arc sec. At z =0.4 for current cosmological parameters (dark 
matter 23%, baryons 4%, dark energy 73%) the distance scale 
would be 5.6 kpc per arc sec. 360 pixels corresponds to 80 kpc. 





R = 360 pixels 



R = 380 pixels 
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FIG. 16: A "neighbors" histogram for clumped galaxies for 
HDF North in the circle with r = 320, 340, 360 and 380 pixels. 



To determine the probability of finding such a sys- 
tematic shift by chance we compute the total number 
of galaxies ( in each of these 300 randomly chosen galaxy 
sets ) for which the circle about a galaxy contains 3 or 
more neighbors. Fig. Eldisplays the results of this study 
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FIG. 17: Histogram of numbers of galaxies in 300 randomly 
chosen sets having 3 or more neighbors in a circle of 360 pixels. 
The red arrow indicates the number of galaxies with 3 or more 
neighbors for the curved set. 

as a histogram: the label of the bins indicates numbers 
of circles with 3 or more neighbors. (The number corre- 
sponding to the "curved" set is indicated by an arrow.) 
This is done for several choices of circle radius. Typi- 
cally, for the optimum radius, which is usually near 340 
to 360 pixels, there will be less than 3 out of 300 sets 
that have as many galaxy-circles with counts equal to or 
greater than the original curved set. In other words, the 
probability of achieving the curved set by chance is equal 
to or smaller than 1%. Since this result holds in both 
fields, and the fields are independent, the probability of 
the observation is less than one in 10*. 

The second method to investigate the clumping is to 
consider a distribution of the distances to the "nearest 
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neighbor" for the galaxies from the "curved" set and com- 
pare it to the random sets (see Fig. I18|l . The horizontal 
axis on Fig. El is the distance to the 'nearest neighbor". 
The Kholmogorov-Smirnov test gives 0.5% probability 
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FIG. 18: A distribution of the distance to the nearest neighbor 
for galaxies in the north HDF. 

that these two distributions are the same. But the situ- 
ation is actually less probable than that, for two further 
reasons. 



^ aligned galaxies R = 350 
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FIG. 19: A "neighbors" histogram for aligned galaxies in a 
circle with r = 360 pixels and the average histogram for ran- 
domly chosen sets in the north HDF. 



Given any distribution of background galaxy orienta- 
tions, the lensing hypothesis would suggest that if the re- 
gions that are populated by dark matter clumps are giv- 
ing rise to the clumping we observe, then equally well, the 
voids should give rise to a population of galaxies that are 
"straighter" , more aligned than normal. In other words, 
we would predict clumping of the aligned galaxies (Fig. 
[T^ ) . We have carried out the same procedure as for the 



"curved" galaxies for "aligned" ones and found that the 
"aligned" galaxies are even more clumped (Fig. I19|l . The 
probability is less than 1 %. 

^ Z distribution for curved galaxies 




FIG. 20: The red (grey) bars indicate the z-distribution of 
"curved" (all) galaxies, respectively . 

^ Z distribution for aligned galaxies 




FIG. 21: The red (grey) bars indicate the z-distribution of 
"aligned" (all) galaxies, respectively.. 

That aligned galaxies are also clumped is not a trivial 
consequence of the fact that curved galaxies are clumped. 
Indeed galaxies midway between curved and aligned are 
not clumped. And though the aligned set is not com- 
pletely independent of the curved set (it represents 1/4 
of the complement of the curved galaxies), we maintain 
that it is independent enough to assert that the probabil- 
ity of finding both clumped by chance is the product of 
the probability of each. With less than a 1% probability 
of curved galaxies being clumped and a 1% probability 
for aligned galaxies being clumped in both the north and 
south fields, the chance probability of all four events is 
less than 1 in a million.'' 

Fig. 1201 shows and compares the z-distribution of 
curved galaxies with the z-distribution of all background 
galaxies. Fig. 1211 compares the z-distribution of aligned 
galaxies with the z-distribution of all background galax- 
ies. Even though the z-distributions for both curved and 
aligned sets are not significantly different from the av- 
erage distribution we have checked the influence of z- 
dependence on the observed clumping. We ran the pro- 
cedures described above for two cases: first, for the ran- 



^ We are not yet claiming the clumping is due to lensing, only that 
there is clumping. 
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FIG. 22: The HDF (North): (a) - triangles (blue dots) are curved galaxies, (b)- boxes (red dots) are "aligned" galaxies and 
black dots are all other background galaxies. 



dom sample restricted to mimic the z-distribution of a 
curved (aligned) set and, second, for z- "blind" random 
sets. The results for the probability of clumping for both 
cases were similar. ® This result was confirmed for both 
fields. We have also carried out an identical process, with 
similar results, for galaxies selected by magnitude rather 
than z-value. As a result of these considerations and ob- 
serving the persistence of clumping under a wide variety 
of circumstances, we are confident that the clumping we 
observe is not occuring by chance. In the next section 
we discuss several possible alternative sources of clump- 
ing. In Fig. \'2'2i a,] we show the spatial location of the 
"curved" galaxies of the north field among all remaining 
background galaxies. In Fig. I22r b) we show the spa- 
tial distribution of "aligned" galaxies of this field with 
all other background galaxies. 



VII. ALTERNATIVE EXPLANATIONS OF 
CLUMPING 

Though we are comfortable in asserting that the 
clumping results cannot occur by chance they may be 
due to reasons other than lensing: 



We have found that when the randomly chosen set is required to 
match the z-distribution of the curved set (in 8 or 10 z-bins) then 
the statistics tests give 1% probability of getting clumping by 
chance compared to 3% probability for the z-independent choice. 



1. Instrumental 

(a) The point-spread function varies over the field 

(b) Pixel derived effects 

2. Derived from the clumping of the background 
galaxies themselves 

(a) Time evolution of background galaxy groups 

(b) Some other group property 

3. Computational 

(a) The coupling of quadrupole and sextupole 
through M2i- 

(b) Galaxy selection and analysis methods 

(c) Image composition (drizzling, overlay, etc.) 

We now address each of these concerns. 

1(a). The point-spread function could presumably 
have a sextupole and quadrupole moment and these could 
be aligned with one another in some smoothly varying 
way across the image. However we feel that this is ruled 
out by the facts that: 

1. the 6 stars in the Hubble north field show no sign 
of having a quadrupole or sextupole strength of the 
required magnitude; 

2. mid-range galaxies (not aligned and not curved) 
could just as well have been clumped under this 
scenario and are not; and 
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TABLE I: Alternative explanations for clumping of curved galaxies. Notations: + => supports hypothesis; OK neutral; 
NO contradicts hypothesis; - =^ contradicts but more data required; NA ^ not applicable. 



Observations 


Lensing 


1(a) PSF 


2(b) Group 
property 


3(a) Mi^ 


3(c) Image 
composition 


1. Clumping of curved 


+ 


OK 


OK 


NO 


- Only along 
boundaries 


2. Clumping of aligned 


+ 


OK 


OK 


NO 


- Only inside 
boundaries 


3. No clumping of 
mid-range galaxies 


+ 


- Could clump 


OK 


OK 


OK 


4. Mixed z content of curved 
clumps 


+ 


OK 




OK 


OK 


5. Curved and clumped concen- 
trated at small a.b. 


+ 


OK 


OK 


NO 

Prefers large a 


OK 


6. Aligned and clumped not con- 
centrated in a,b. 


+ 


OK 


OK 


OK 


OK 


7. Direction to scattering center 
varies 


+ 


- 


OK 


OK 


NO 

Pattern expected 


8. Curved galaxies next to not- 
curved galEixies 


+ 


NO 


OK 


OK 


NO 

Pattern expected 


9. Stars are round 


OK 


Not enough stars 


OK 


OK 


Not enough stars 


10. Cluster CL0024 


+ 


NA 


NA 


NA 


NA 


11. Deduced mass magnitudes 
have reasonable values 


+ 


NA 


NA 


NA 


NA 



3. the deduced directions to the scattering center are 
erratic. 

The logic behind 2) is that if the point-spread func- 
tion is causing this effect, its sextupolc^ moment and 
quadrupole moment are continuously varying across the 
image so as to induce the observed distortion in the im- 
ages. There could be isolated areas of the sky were this 
effect is mid-range (not necessarily aligned or curved). 

1 (b) . We have carried out studies to see if unexpected 
moments are generated by dividing an image into pixels. 
For example, we took a known bi-Gaussian distribution 
and, varying the centroid, projected it onto a pixel grid. 
The falsely induced sextupole moment had strengths less 
than 10^"''. In general, there is no reason we have dis- 
cerned for which pixelation effects lead to spatial corre- 
lations, since pixels (modulo pixel defects) are uniform 
across the field. 

2(a). Since the galaxies at a slice in z are known to 
be clumped, then if there is some age-dependent change 
in galaxy shapes, a shape selection criteria could be see- 
ing an age-biased sample which could then be clumped. 
Since our "curved" galaxy set has the same i;-distribution 
as all galaxies, the premise would appear to be false, ie 
galaxy curvature as we are quantifying it does not appear 
to evolve with time. 

2(b). The premise is that galaxy groups possess some 
property (other than age) that effects the shape of galax- 
ies. Perhaps there was an explosive event or set of events 
that extended across the group (80 kpc radius) that ef- 
fected the birth process of the background galaxies. If 
true, it would be an unexpected and interesting result 
in itself. We don't think we can rule this out yet, but 



with more statistics an argument to discount this sce- 
nario could be based on the fact that the curved and 
clumped galaxies in any given clump have a variety of 
z- values. 



3(a). If is non-zero, a AM20 can lead to a corre- 
lated quadrupole and sextupole moment. A clumping of 
"curved" galaxies could result from a spatial dependence 
of the orientation of with respect to AM20. A non- 
lensing origin of such a spatial dependence is an alterna- 
tive source of clumping, but would require its own expla- 
nation, presumably one of the other alternatives among 
the items in our list. 



3(b). Galaxy selection would appear to be blind to 
position in space and incapable of leading to the spatial 
clumping of curved or aligned galaxies 



3(c). This represents our most serious concern. We 
suspect that we are processing this data in a way which 
was not anticipated by the creators of the image construc- 
tion process, and we are concerned, for example, that on 
the boundaries of overlays, distortions could be present 
that would give the images a false shape. There are not 
enough stars in the field to rule this out using star im- 
ages. Our best response to this concern is to note that 
the regions where curved galaxies clump do not appear 
to have any particular identifiable pattern, i.e. they do 
not appear to coincide with sub-field boundaries. 
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a-b plof for curved and all background galaxies 
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FIG. 23: The HDF (North): triangles correspond to "curved" 
both plots are all other background galaxies. 



VIII. THE a b MOMENT-MAGNITUDE 
PLANE 

Fig. I^ni (left) is what we call an ^a, bj plot for the 

sample of all galaxies (dots) and the set of curved galax- 
ies (triangles). Each point on the plot corresponds to a 
single galaxy. The horizontal axis is the magnitude of 
the deduced sextupole coefficient b, and the vertical axis 
is the magnitude of the deduced quadrupole coefficient a. 
Note that galaxies having both large a and large b which 
occur in the set of all galaxies are noticeably absent in the 

set of curved galaxies. In the ^a, b^ plot for all galaxies 

40% have both a and 6-values below the mean, whereas 

in the ^d, plot for curved and clumped galaxies, 70% 

have both a and b values below that same mean. 

This trend was predicted for lensing, as a consequence 
of considering the vector addition of the original moment 
and the induced moment. When the original moment 
is larger than the induced moment, the resulting vector 
tends to be aligned in the direction of the original mo- 
ment, and since this angle is then divided by 2 in the case 
of the quadrupole moment, and by 3 in the case of the 
sextupole moment, the alignment with the original mo- 
ment is much better than would be normally expected. 
On the other hand when the induced moment is larger 
than the original moment the argument works the same 
way to deduce that the alignment with the induced mo- 
ment is much better than expected. 

In other words, the induced curving is expected to 
be seen if the background moments are as small as, or 
smaller than the induced moment. So the population on 
the a — b plot for curved and clumped galaxies is a strong 
indication of the strength of the induced moments. Of 
course, any other "add-on" effect will have the character- 
istic that it will be more noticeable for small normal or 
original moments, but the vector addition law with the 
29 and 39 dependence has a remarkably sharp behavior. 

"Aligned" galaxies, on the other hand, are presumed to 
represent galaxies that have not been altered by lensing 



a-b plof for aligned and all background galaxies 
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(left plot) and "aligned" (right plot) galaxies and black dots on 



into a mid-range or curved shape. Thus their distribution 
on the a-b plane should be very similar to all galaxies. 
This is indeed the case, as can be seen in Fig. (right 
plot). 

IX. ESTIMATE OF GROUP MASS 

The results of this and the following section are more 
speculative. We assume that the observed clumping is 
indeed due to lensing and attempt to deduce the proper- 
ties of the lensing mass distribution. This is premature 
because we have not done the modeling to determine sys- 
tematics and because the sample of clumped and curved 
galaxies is small (total of 110, both fields). On the other 
hand, we feel it is possible and appropriate to make order 
of magnitude estimates. 

An estimate for the mass Mx of the group could be 
written 

Mx^pAx = ^Ax = ^(^) Ax, (26) 

where M is the mass of the constituent clumps. Ax the 
area of the group, and dx the typical separation distance 
of the clumps within the group, rg is chosen so that in- 
sertion of M/vq into the formula for a produces a typical 
induced moment size for events that change the popula- 
tion of curved galaxies. Pa = irrQ/d^ can be interpreted 
as the probability that any particular light path receive a 
moment change of the required magnitude. This can be 
found by obtaining an estimate for the probability that 
the light path of any particular galaxy passing through 
the group receives noticeable induced moments. 

To obtain this probability we have attempted to iden- 
tify the clumped parts of the field (which is justified since 
we know it is indeed clumped) and determine the frac- 
tion of the galaxy light streams which become curved 
when penetrating these parts. We have done this using 
a friends-of-friends algorithm, requiring that i) there are 
at least 4 curved members in a group, and ii) that to be 
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included in a group a curved galaxy must lie within a 
certain distance (about 350 pixels) of other members of 
the group. Fig. |^is a combination of fig. I22f a') and fig. 
|22Ib) with the "curved" and aligned galaxies identified 
as belonging to clumps connected by straight lines. We 
note that we have groups with area ranging from 6 10"' 
pixels to 1.5 10^ pixels. Within the groups about 40% 
of the galaxies are curved. Outside of the groups less 
than 10% of the galaxies are curved. So the probability 
of being transformed to curved is estimated at 30%. 

The transformation to curved depends on the suitabil- 
ity of the background galaxy (for example its moments 
are small enough to easily alter) times the probability 
that the kick is large enough. We assume that these are 
roughly equal, and hence get a probability of Pa ~ 0.6 
that the kick is large enough. 

This yields the estimate 



MxG = 



O-typ DtS 

47r DtlDls 



PaA 



X 



(27) 



this degeneracy because the expression for the induced 
quadrupole moment of equation ^ determines only the 
ratio M/rQ. Estimating the clump mass requires addi- 
tional information. 

Additional information can come from an analysis of 
the relative role played by the quadrupole and sextupole 
moments, according to the following considerations. It 
follows from equation © that the strength of the induced 
sextupole moment occurring along with the induced 
quadrupole moment atyp would be btyp = atyp rc/ro. 
Assuming the lensing plane is at z = 0.4, re ~ 0.5 kpc 
and since Ohjoa ~ 1/5, and as we will argue below that 
ajh K, a/h, then both sextupole and quadrupole moments 
will be equally observable for an tq « 2.5 kpc. Since we 
have estimated above that Pa = Tfr^/dx ~ 0.6, for this 
ro we have d ~ 6 kpc and N = Ax/d\ w 500 for the 
typical groups. Thus for htyp = 1/5 atyp the estimate for 
the clumping mass would be M m 10^" Mq. Following 
this line of reasoning we have constructed the following 
table. 



w (0.5 to 2) atyp[]yc]. 

atyp has to be large enough to change moment align- 
ments, so it must be comparable to but probably less 
than a a- Also we must account for the systematic weak- 
ening of moments by the point-spread function, which 
(see below) we would estimate at a w 2d. For atyp ~ 0.2, 
MxG w (0.1 to 0.4) pc^Mx = (2 to 8 lO^^) 

Moments a cannot be larger than aa because the dis- 
tribution of moment strengths for curved galaxies is ac- 
tually smaller than that of all galaxies. The realignment 
of moments can be thought of as a change of direction of 
the quadrupole moment or a change of direction of the 
sextupole moment, or a combination of both. 

A distinct possibility is that there is no induced sex- 
tupole; only the quadrupole moment direction is changed 
leaving a sextupole behind resulting in a "curved" galaxy. 
One might think that higher order moments get moved 
along with the quadrupole because the map is linear. For 
example a bi-Gaussian, which has its octupole moment 
aligned with the ellipticity, is mapped into a bi-Gaussian 
which will also have its octupole moment aligned with 
its ellipticity. However, for sextupole moments this ques- 
tion is answered quantitatively by equation H12I) which 
includes a term showing that with 6 = 0, the observed 
sextupole moment change will depend on the magnitude 
of Mji. 



TABLE 11: Possible lensing parameter range. 



Lensing mass /A/© 




10^" 


10^^ 


10^^ 


Number of clumps in 


5000 


500 


50 


5 


largest groups 










Typical impact parameter 
ro (kpc) 


0.8 


2.5 


8 


25 


b/a ratio 


0.6 


0.2 


0.06 


0.02 


dx (kpc) 


2 


6 


20 


60 



From a lack of a population of localized large b values, 
masses below 10^ are ruled out. Under more careful in- 
spection and better statistics the range from a few 10^ 
through 10^^ should have observable consequences for b. 
The presence of lensing masses greater than 10^^ would 
be hard to numerically constrain from looking at sex- 
tupole distributions, except to know the mass is larger 
than the observable cut-off. However, there should be 
observable consequences for the quadrupole distribution, 
since there would be a population of larger induced mo- 
ments. 

In other words, we believe that a careful study of the 
distributions of a and b together with modeling should be 
able to provide sufficient statistical information to deduce 
the distribution of lensing masses. 

We end this section with a discussion of systematic and 
random error. Important error and noise sources in any 
individual measurement include: 



X. ESTIMATE OF CLUMP MASS 

The typical groups in the field could presumably be 
composed of 5 clumps of mass lO^^M©, 50 clumps of 
mass 10"Mq, 500 clumps of mass IO^'^^Mq, 5000 clumps 
of mass lO^^M©, or some fractional combination of these. 
A single mass at 5 10"'^^ solar masses is ruled out by the er- 
ratic directions to the scattering center. We are left with 



1. the systematic effects of the point-spread function 
and thresholding, 

2. random background galaxy moments, 

3. extraction of moments from the image, 

4. pixelation of images, and 

5. photon counting noise. 
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FIG. 24: The HDF North red boxes are "ahgned" galaxies, blue triangles are curved galaxies and black dots are all other 
background galaxies. 



(1) The largest effect arises from thresholding and the 
PSF. For our galaxy cores, we estimate MH « Mjj^, 
whence Ki ^M^^, implying a ^ 2a. For Gaussian 
shapes, M22 = 4A/^j^, hence an estimate for the denom- 
inator of the expression for b is approximately 10 cr* ^ 

2.5 M22 and the numerator has an Af ^ •\/2M^, 

yielding the estimate b ^ 2b. These are large adjust- 
ments that need careful attention, but the gross indica- 
tion is that the a/b ratio will suffer a smaller adjustment 
than either numerator or denominator and led to the use 
above of an estimate a/b w a/b. 

(2) The next largest uncertainty in the a/b ratio comes 
from the existence of non-zero background moments. 
It is a happy circumstance that there is clumping of 
"aligned galaxies" as well as "curved" galaxies, because 
the clumps of aligned galaxies are presumed to indicate 
a lack of lensing within their domains. Hence the dis- 



tribution of galaxies behind voids will give an important 
base to which the galaxies behind lensing clumps can be 
compared. With larger statistics, one could expect to 
extract interesting details by studying these differences. 
Also it will be invaluable to have a precise knowledge 
of the point-spread function, including all its moments, 
throughout the field. 

(3) We feel that our method to determine moments, 
validated by its ability to reveal important correlations, 
can be substantially improved. This is intimately related 
to item (1). 

(4) In pixelation studies, we were surprised to find that 
the change of the centroid distribution was unable to pro- 
duce apparent sextupole moments above the 10~^ level. 
Square pixels do give rise to spurious octupole moments, 
and some care is required in that case. 

(5) Our thresholds are typically set at approximately 
175 photons per pixel. The core peaks are the order of 
600 photons counts per pixel. We have not studied the 
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effects of photon counting noise, however we would not 
expect its obvious random nature to affect our conclu- 
sions, and it should be substantially smaller than item 
(2). 



XI. DETERMINATION OF CLUMP RADII 



The probability of penetrating a clump would go like 



Probability = ~ 



M ri 



M 



X 



M 



'0 "x 



M T:r\' 



(28) 



where is the lensing mass (clump mass) radius. The 
last ratio in this equation is fixed, so the probability 
depends only on the first ratio (the inverse of the pro- 
jected density within a clump) which could be slowly 
varying. While one expects the 3-dimcnsion density of 
small clumps to be larger than large clumps, the pro- 
jected densities can be similar. So the probability of 
penetrating a clump could depend very weakly on the 
composition of the group. 



relative 
strength 



Moment strengths 

iquadrupole 




FIG. 25: Strength of the quadrupole, sextupole, octupole, 
and cardioid-like kicks within a Gaussian clump. The hori- 
zontal axis is the radius divided by the RMS lens radius. 

Penetration within a clump mass distribution will af- 
fect the strength of the induced moments and should 
leave a signature on the sextupole and quadrupole mo- 
ment distributions. This effect is larger than one might 
initially expect because the sextupole moment within a 
clump is strongly suppressed and only reaches its asymp- 
totic form at r = Sr^ . Figl^S] illustrates the situation 
fully for a Gaussian mass distribution. If > vq [tq 
is the background galaxy light-stream footprint radius) 
then the interior can be well probed, and the quadrupole 
to sextupole ratio will remain less than 

h 1 ra 1 
- < -— < - • 

a 3 tl 3 

If tl < ''g/S, then as tq re the sextupole to 
quadrupole ratio can approach 1. In other words, for 



both a and 6, one must compare the difference distri- 
butions between voids and clumps. The behavior of the 
tail at large moments will have a different distinguishable 
behavior depending on the radii of the lensing masses. 

The interior may also be probed with the d coefficient, 
which is zero except within distributions. Its strength 
distribution is also shown in Fig. 1251 Over its limited 
range (the maximum is at rg = tl) its value is surpris- 
ingly large. The expression for finding d from moments is 
given by equation H15|l . Evidence suggesting that galaxy 
light paths are penetrating lensing mass distributions is 
present in the angular asymmetry evident in Fig. 1141 
The background distribution of d would be expected to 
be symmetric. 



XII. SUMMARY 

We visually examined faint images selected by the SEx- 
tractor software from the Hubble deep fields using an un- 
usually high threshold. After filtering images with two 
or more maxima, we measured sextupole and quadrupole 
moments. The "curved" galaxies we sought were iden- 
tified as those whose sextupole moment was oriented so 
that one of its minima was aligned within a few degrees of 
a quadrupole minimum. Wc then looked for and found an 
improbably large spatial clumping in each Hubble deep 
field of both curved and aligned galaxies. The probability 
of each is the order of 1%. 

Our motivation and preferred hypothesis is that these 
galaxies were lensed by close collisions with, for example, 
10"'^° solar mass clumps that reside within a half dozen 
groups of mass a few times 10^^ solar mass in each field. 
The projected spacing of such clumps within these groups 
would be about ^.5 kpc. The rms radii of the footprint 
in the lensing plane of the observed background galaxies 
is about 0.4 kpc, a factor of 10 smaller than this spacing. 
The cores of these galaxies could well be the order of 2 
kpc, so there is a hope that the effects of light paths trav- 
eling through clump interiors may be observable through 
the d coefficient introduced above. 

We have carried out other tests that support the hy- 
pothesis that the observed clumping comes from lens- 
ing: 1) aligned galaxies were predicted and found to 
be clumped as strongly as curved galaxies, 2) galax- 
ies halfway between aligned and curved were predicted 
and found to have no clumping, and 3) correlations were 
predicted to be more readily detected for galaxies with 
smaller moments. 

We have constructed alternate hypotheses for our ob- 
servations based on instrumental effects, computational 
effects, or other physical phenomena. The tests we have 
constructed and their implications for each alternative 
hypothesis were discussed and the results summarized in 
Table I. These alternatives can be eliminated (or con- 
firmed) with more data. 

Finally, the numbers that we are seeing have very in- 
teresting and plausible magnitudes, they even may be 
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expected in the standard structure formation scenarios 
There are four important future studies: 

1. Understand the image construction process of the 
deep-field images to insure that the clumping prop- 
erty is not an artifact of this process; 

2. Simulations of small- impact parameter lensing to 
quantify the relationship between moment distri- 
butions and mass distributions and to discern the 
fraction of galaxies that are lensed but not seen; 

3. Improvement of the image-analysis process, includ- 
ing variable thresholds and point-spread function 
removal; 

4. Deep field studies with the additional pictures from 
the ACS camera on Hubble that will be available 
in the coming months and years. 

In truth our observations should not be considered to 
be a subtle effect. Our samples from each Hubble deep 
field contains only the order of 350 background galax- 
ies satisfying our selection criteria. About 75 of these 
are "curved" , and of these about 60 appear to reside in 
clumps. In other words, some 16% of background galax- 
ies would be experiencing observable "small impact pa- 
rameter" scattering from dark matter clumps. 

These results are possible because i) the background 
sextupole strengths are a factor of 4 smaller than the 



quadrupole moments, ii) the sextupole and quadrupole 
moment orientations arising from lensing are correlated, 
and iii) one can look for clumping on the sky. 

The Hubble deep fields are less than 2 min by 2 min, 
about 10"'^ sq. deg. Projects in planning stages (eg. 
SNAP) have a weak lensing program of between 300 and 
1000 sq. degrees with resolution comparable to the Hub- 
ble deep field observations. 

If our conjecture were to be true, weak sextupole lens- 
ing could provide valuable insight into mass structure at 
length scales 100 times smaller than weak quadrupole 
lensing, at smaller scales than the Lyman alpha forest 
results 01 . 
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